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^f) . The models of large extra compact dimensions, as suggested by Arkani-Hamed, 

■ Dimopoulos and Dvali, predict exciting phenomenological consequences with gravita- 

Q^' tional interactions becoming strong at the TeV scale. Such theories can be tested at the 

Qh' existing and future colliders. In this paper, we study the contribution of virtual Kaluza- 

r~| , Klein excitations (both spin-0 and 2) in the process e'*"e — > HHZ and e~^e — > HHj 

^> . at future linear collider (NLC). We find that the virtual exchange KK gravitons can 

^ ■ modify the cross-section a{e~^e'~ — > HHZ) significantly from its Standard Model value 

' and will allow the effective string scale to be probed up to 6.6 TeV. The second process 
is absent at the tree level in the standard model, and can therefore be used to put 
limits on the effective string scale of 7.4 TeV. 
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1 Introduction 



The concept of large extra dimensions and TeV scale gravity introduced by Arkhani-Hamed, 
Dimopoulos and Dvali, better known as ADD model [1] has attracted a lot of attention. 
In this scenario, the total space-time has D — A + n dimensions, where gravity lives in n 
additional large spatial dimensions of size R. The standard model particles are confined to 
the usual 3 + 1-dimension. 

The effective AD Planck scale, Mpi ~ 2.4 x 10^^ GeV , is related to the fundamental 
Planck scale Ms in (4 + n) dimension by 

~ i?"M^+2 (1) 

Thus, for the large extra-dimensions, it is possible to have a fundamental scale Ms as 
low as a TeV [1], solving the gauge hierarchy problem of the standard model. With Ms ~ 
TeV, for n = 1 the value of R (~ 10^^ millimeter ) is ruled out by gravitational experiments. 
On the other hand, n > 2, we get R < mm, a range that is still allowed by gravitational 
experiments. 

Each graviton couples to the standard model matter and gauge particles through energy- 
momentum tensor and its trace, with the strength suppressed by powers of the 4-dimensional 
Planck scale, Mpi. However, from the 4-dimensional point of view, the massless gravitons 
propagating in the (4 + n)-dimensional bulk are seen as massive towers of Kaluza-Klein (KK) 
modes of excitations with spin-0, spin-1 (which decouples ) and spin-2. The mass spectrum 
of these KK modes can be treated as a continuum, since their mass splitting ~ 1/Ris about 
10~^ eV for n = 2 and for n = 6 it is of the order of a MeV. After summing over these 
KK modes one gets enhancement of KK graviton coupling to standard model matter and 
gauge particles to powers of {^/s/Ms) [1], where ^/s is the available energy for the process. 
The Feynman rules for this theory have been developed considering a linearized theory of 
gravity in the bulk in Ref. [2, 3]. These new interactions can give rise to several interesting 
phenomenological consequences testable at present and future coUiders [4]. The effect of 
these new interactions can be observed either through production of real KK modes, or 
through the exchange of virtual KK modes in various processes [4, 5, 6]. 

In this paper we will consider the process e'^e~ — > HHZ and e'^e~ — > HH^ to study the 
effect of low-scale gravity at proposed linear colliders with center of mass energy 500 GeV 
and beyond. However, one should keep in mind that these processes will not serve as the 
dominant discovery channel for KK gravitons, since there are more direct processes for its 
discovery [4]. Once the discovery of KK gravitons is established, then the next phase will be 



to look for their effect on some more complicated processes, like the ones discussed here, to 
re-confirm their discovery. 

In the standard model, HHZ production has been studied in the context of the determi- 
nation of the Higgs self couplings [8] . The other channel, HH^, is a unique process mediated 
by virtual gravitons. This process is absent at the tree level in the standard model, there- 
fore, the observation of such a process can lead to a possible indication of low scale gravity. 
Through out this paper we assume that the Higgs boson will already have been discovered 
and its mass determined. 

The rest of the paper is organized as follows: In Section 2, we study the production 
of HHZ and discuss the effect of KK gravitons. In Section 3, we will discuss how virtual 
exchange of gravitons give rise to HH'-f final state and discuss the consequences. Section 4 
is reserved for the overall discussions and conclusions. In Appendix, we present the graviton 
exchange amplitudes for e"'"e~ —>■ HHZ process. 

2 e+e- ^ HHZ 

In this section we study the additional contribution of virtual exchange of spin 2 and gravi- 
ton modes in models of large extra dimensions to the Higgs pair production in association 
with Z. In the standard model, the process e+e" HHZ has been studied as a means of 
determining HHH couplings. The standard model amplitudes and polarized cross-sections 
are given in the Ref . [8] . In our analysis, we shall use unpolarized beams. Additional diagrams 
arising from the exchange of both spin-2 and spin-0 KK gravitons are shown in Figure 1, 
where KK stands for both the spin-2 and spin-0 modes. The amplitudes using the Feynman 
rules ^ given in Ref. [3] are presented in the Appendix. 

After adding coherently these new amplitudes to the standard model ones, we get the 
total amplitude. The new amplitude now depends on Higgs boson mass(Mij) , center-of-mass 
energy (y^), fundamental Planck scale (M5), and finally on the number of extra dimensions 
(n). Through out our analysis we have fixed the standard model Higgs boson mass at a 
representative value Mh — 120 GeV. For center of mass energy, we consider three possible 
values of ^/s — 0.5 TeV, 1 TeV and 3 TeV at which the future e^e" colliders are expected 
to operate [5, 6, 7]. 

In Figure 2, we represent the total cross-section for Mh = 120 GeV and n = 3 as a 

^In our analysis we take ^ = 1, the de Bonder gauge. 



function of machine energy for three different values of Ms — 3.5,4.5 and 5.5 TeV. The 
dashed hne represents the standard model prediction and this exhibits the expected fall-off 
behavior with energy. For large values of Ms, this behavior is preserved, since the graviton 
contribution then becomes small. However, when Ms is smaller, the cross-section show a 
marked increase with energy, showing the well-known feature of gravitational interaction. 
It is obvious from this Figure, that at energies around 3 TeV, the graviton contribution 
is enormous, if Ms is small ( 3.5 TeV); however, a discernible difference exists even when 
Ms = 5.5 TeV. Thus we can expect larger effects, or in other words stronger bounds on Ms 
as the machine energy increases. 

In Figure 3, we represent the total cross-section as a function of Ms for different values of 
extra dimensions n = 3, 4, 5, 6 for two value of machine energies, ^/s = 500 GeV and 1 TeV. 
For comparison the standard model cross-sections are given by the dashed lines, which are 
independent of Ms- From this Figure, we can sec that for a given machine energy and Ms, 
the graviton contribution decreases with the increase in extra-dimension (n). For a given n, 
the graviton contribution asymptotically tends to the standard model value with increase of 
Ms as it should. 

In Figure 4, we represent the angular distribution of the Higgs boson. The standard model 
distribution (dashed line) shows typical sin^ 9 behavior, whereas, the ADD contribution 
(solid line), which has in addition spin-2 and spin-0 exchange, can be distinguished from the 
standard model very easily. From this distribution, it is evident that larger deviation from 
the standard model is expected in the forward and backward regions of the detector, with 
0.5 <| cos^ |< 0.8. We assume that the efficiency of the detectors in that region will be 
good enough to observe such events. In the central region ( | cos^ |< 0.2 ) of the detector, 
the standard model cross-section shows a peak over the ADD contribution. 

In Figure 5, we represent the fractional deviation R — ^'^^^^^^'^^'^^ as a function of Ms for 
two values of center of energies, 500 GeV (dashed line), and 1 TeV ( sohd hne ), where, a (Ms) 
is obtained from coherent sum of ADD and standard model amplitudes. The different choices 
of extra dimensions (n) are shown along each curve. It is very clear from this Figure, that 
higher the machine energy, larger is the deviation from the standard model. We consider that 
more than 10% deviation (i? = 0.1) from the standard model prediction with an integrated 
luminosity C of 2ab~^ will be able to provide reasonable evidence of low-scale gravity. This 
argument is based on the result obtained in Ref.[9]. In this paper, the authors have studied 
the standard model ZHH production at future e~^e~ linear collider. They have shown that 
at y/s = 500 GeV, the e+e" — > ZHH cross-section can be determined with a 10% relative 



error assuming an integrated luminosity of 2ab~^. In our analysis we use this error (10%) 
on the cross-section measurement to determine 5a discovery limit on Ms with the same 
integrated luminosity. 

For higher center-of-mass energies(y^ = 1 TeV and 3 TeV), we also assume for simphcity 
the same relative error in the ZHH cross-section measurement with the same integrated 
luminosity. We think as a first exercise of this kind, the assumption considered here is 
sufficient. Before we display our limits on the scale Ms, we would also hke to stress that the 
limits on the string scale Ms obtained here are merely indicative. These limits may change 
when one considers detector simulations with full estimations of signal and backgrounds. 
As we have already mentioned, this process should not be considered as the main channel 
to probe the low scale quantum gravity. Since, by the time the next future e^e~ collider 
starts operating and reaches such a high integrated luminosity, the low scale quantum gravity 
model would have been discovered if true. 

In the Table 1, we display the 5a discovery limits on Ms as a function of e^e~ center 
of mass energy number of extra-dimensions (n) and integrated luminosity. From this 

Table, we can appreciate the following features: 

• The 5(7 limit on Ms is weekly dependent on number of extra dimensions (n). 

• The 5(7 limit on Ms obtained for 500 GeV machine is not very promising. In this case 
the strongest hmit is less than 1 TeV for n — ?>. 

• The situation become better once we go to a higher center of mass energy. The strongest 
limit on Ms is 6.6 TeV obtained for n = 3 and at ^/s = 3 TeV. 

3 e+e- HH-f 

For the completeness of our study, in this section we discuss the production of HH^ at e+e~ 
collider through the virtual exchange of spin-2 and spin-0 KK gravitons. The scalar pair 
production at e+e~ collider in the models of large extra-dimensions has been studied [10]. 
The Feynman diagrams for e+e~ — > HH^ are similar to Figure 1, with Z replaced by 7. In 
the case of exchange of spin-2 gravitons all the four diagrams (a — d) contribute, while for 
spin-0 KK mode exchange, only diagrams (a), (6) and {d) contribute. These three diagrams 
form a gauge invariant set, while the diagram (c) vanishes due to Dirac equation and gauge 
invariance. In this process, unlike Higgs boson pair production, spin-0 mode of KK graviton 
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Table 1: 5a limit on Ms assuming an integrated luminosity of 2 ab ^ for different center of 
mass energies {-s/s). n is the number of extra dimensions. 



also contribute. The photon will serve as an additional trigger. We impose the following 
identification criteria for the photon [1 1] : 

> 10 GeV (2) 

\cose^\ < 0.98 (3) 

where, 9^ defines the photon angle with the e" beam direction. This angular cut on the 
photon removes the coUinear divergence. These 'acceptance' cuts are more-or-less the basic 
ones. Though further selection cuts will become appropriate when a more detailed analysis 
is done, it suffices for our analysis, which is a preliminary study. 

Unlike the previous case, here the whole contribution comes from the exchange of virtual 
gravitons. As before, we fix the Higgs boson mass at the value used in the last section as 
also the machine energies. 

In Figure 6, we represent the total cross-section for Mh — 120 GeV and n = 3 as a 
function of machine energy for three different values of Ms — 3.5,4.5 and 5.5 TeV. For a 
given Ms, the cross-section shows a marked increase with energy, depicting the well-known 
feature of gravitational interaction. 

It is obvious from this Figure, that at energies around 3 TeV, the graviton contribution 
is enormous, if the Ms is small ( 3.5 TeV); however, a discernible difference exists even when 



Ms — 5.5 TeV. Thus we can expect larger effects, or in otfier words stronger bounds on 
Ms as tfie macliine energy increases. In Figure 7, we represent tfie angular distribution of 
the Higgs boson. This particular shape of the distribution is specific to virtual exchange of 
gravitons in the process. 

In Figure 8, we represent the total production cross-section as a function of Ms for 
different values of extra-dimensions n = 3, 4, 5, 6 for three values of center-of-mass energies, 
^/s = 500 GeV (dotted hnes), 1 TeV (dashed lines ) and 3 TeV ( sohd lines). For a given s/s 
and Msi the graviton contribution decreases very slowly with increase of extra-dimension n. 

The main possible source of background to this process is e^e~ ZH'-f with Z decays to 
bb looking like a Higgs decay. Using the CompHEP [12] we have computed o'^/(= o'{e~^e~ 
ZH'-f) X Br{Z bb)) x Br{H bb)) for e^e~ ZH'-f process for three values of = 0.5 
TeV, 1 TeV and 3 TeV to get a feeling for this background. The numbers are 0.50 fb, 0.15 
fb and 0.023 fb for = 0.5 TcV, 1 TeV and 3 TcV respectively. There are two ways one 
can eliminate this background. Firstly, if the Higgs mass is heavier than Mz (which is in 
fact true from LEP 11 data) [13], the invariant mass distribution of particles produced from 
Higgs decay will be able to reduce the background. Secondly, it should also be possible to 
reduce the background further by studying the angular distribution associated with these 
processes. The loop-induced standard model or MSSM Higgs pair production cross-section 
for hght Higgs mass is of the order 0.1 — 0.2 fb at = 500 GeV [14]. This cross-section 
will be suppressed by order ccem, if an additional photon is produced. When folded with 
H ^ bb branching ratio, the cross-section becomes 0(10"^) fb. This cross-section will be 
further suppressed by the 3 body phase-space. 

To estimate the possible 5a discovery limit on the scale Ms with an integrated luminosity 
of 2ab~^, we only consider the ZH^ as the main source of standard model background. We 
multiply the standard model cross-section by the luminosity to get the predicted number of 
events. We then estimate the errors assuming that the statistical errors are Gaussian and 
that there are no systematic errors. This certainly makes our estimates of the discovery limits 
over-optimistic ( especially for for ^/s — 3 TeV, where number of standard model events are 
small ). In any case, before more detailed studies of the standard model backgrounds and 
the detector design, any estimate of errors must be considered a crude estimate. In Table 2 
we display such discovery limits on Ms for different machine energies and extra-dimensions. 
As expected the strongest limit appears for ^/s = 3 TeV and for n = 3. 
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Table 2: 5a discovery limit on Ms assuming an integrated luminosity of 2 ab ^ for different 
center-of-mass energies {-s/s). n is the number of extra dimensions. 

4 Conclusions 

In this paper we have studied the implications of KK graviton contribution to the process 
e+e~ — > HHZ, which has been studied in the standard model. The spin-2 and spin-0 
mode of KK gravitons contribute to this process substantially. However, the existence of 
low-energy quantum gravity may be discovered through more direct channels as studied by 
several groups. Nevertheless, this process will be an independent confirmation for such a 
discovery. The 5(j discovery limits on the string scale Ms is obtained assuming an integrated 
luminosity of 2 ab~^. We have derived this limit considering 10% error on the HHZ cross- 
section measurement at ^/s = 500 GeV at £ = 2 ab~^. In reality for higher center-of-mass 
energies, the above error may change. However, in our analysis we have assumed it to be 
remain the same. The bounds obtained here are only suggestive and may change when 
one take into account proper detector simulation, inclusing H and Z decays and imposing 
selection cuts on the final state particles and the full background calculation. 

Another interesting behavior is shown by the Higgs angular distributions. In the standard 
model, the shape behaves like ~ sin^ 9^ whereas it is completely different in the presence of 
KK gravitons. Furthermore, we have shown that in the forward and backward regions of the 
beam pipe, one would expect larger deviation from the standard model prediction. Hence, 



careful study of Higgs angular distribution may provide evidence of new physics beyond the 
standard model. 

To complete our analysis we have also studied the e+e~ — > HH^ process in this model 
which is absent at the tree level in the standard model. We have discussed how this mode 
can be distinguished from the possible backgrounds. We have obtained the Scr discovery 
hmit on Mg assuming an integrated luminosity of 2 ab~^. 



5 Appendix 

Here we will write down the graviton contribution ( both spin 2 and ) amplitudes for the 
process :e+e~ — > HHZ. In our analysis we have chosen ^ = 1, the so called de Donder gauge. 
The Feynman rules used to obtain these matrix elements have been taken from Ref . [3] . The 
relation between the fundamental scale Ms and the size R of the n extra dimensions is given 

by [3] 

K^R" = 167r(47r)"/2r(n/2)M5 (4) 

where, k = x/IOttGat. Gat is the four dimensional Newton constant. The amplitudes a — d 
and e — h correspond to the exchange of a virtual graviton of spin 2 and respectively. 

Ma = -tJ-Gjt r^X^p[v{k2)C^.{h-^ih''{9v+9aj5)u{h)]e;{pi) (5) 

ZCw (fci — Pi) 

Mt = i^GjT 7^X^0[Hk2h%9v+9al5){h-^>iKMki)]e;ipi) (6) 



2Cw {k2-pif 



= C^^ {s - Mlf '^^'^^^''" ~ ^)X.Mk2)l,{9v+9al^)u{k,)\e^*{pi) (7) 



w 



M, = ^(_iL_)^(p2.p3 + 2Mi)[z;(fc2)7"(fl.+5a75)«(fci)]e*(pi) (9) 

Mf = Me (10) 

= - (^S^^ {-J+lt- Ml) ^P'-P' + '^^)^^^ Hk2h.{9. + 9al,)u{k,)] {g'^^ - ^)(11) 

^'^ = f^(z;^)^(f2.P3 + 2Mi)[^;(fc2)7"(5.+5a75)«(fci)]e:(pi) (12) 

where, q — {ki + k2), ki and k2 are two incoming four momenta. Cw — cos^vk, dw the 
Weinberg angle, and Qa are the vector and axial vector couplings of Z to electron. 

a = (13) 



XaP = {M]jriap - CapSvPiPs) (14) 

2|3^^ = {r]p^Ml + q'f^q'^ - pi^q'f^) (15) 

Cf^u = bf^ih -pi- k2)u + (yu ^ i^)] - '^Vnyih -pi- h) (16) 

= C^.(pi^-Pi) (17) 

■^pai^A = [{Ml - q.pi)Cp„^^x + Dp„^^x{.(l,-Pi) + Epa^^x{.q,-Pi)\ (18) 

pM'^a/? ^ ^t^a^up ^ _ ^_^^/xz.^a/3 Donder gauge) (19) 



where, si = (^2 + ^3)^, P2 and ps are four momenta of outgoing Higgs particles. The function 
D[si) counts for the virtual KK state exchanges. Complete expression of D{si) can be 
found in Ref. [3]. q' — and rj^'' is the Minkowski metric tensor. The definitions of 
Cpa,ux, Dpcr,^x{q, —pi) and Epcr^i,x{q, —pi) can be found in [3]. The expression for pf^'"^!^ is 
taken from Ref. [15]. The amplitudes for the process e^e' — > HH^ can be obtained from 
Equations (5-12), by putting Mz = 0, = 1 and Qa = 0. For spin-2 graviton mode exchange, 
all the four amplitudes in Equations (5-8) will contribute, while for spin-0 mode. Equations 
(9,10) and Equation (12) will contribute. 
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Figure 2: Variation of the cross-section (fb) for HHZ production with machine energy. The 
sohd curves correspond to the ADD predictions for Ms — 3.5, 4.5 and 5.5 TeV and n — ?>, 
while the dashed hne represents only the standard model contributions. The Higgs mass is 
set to 120 GeV. 




Figure 3: Variation of the cross-section (fb) for HHZ production with Ms for = 0.5 and 
1 TeV are shown by the sohd hnes. The values of extra dimensions (n) are given along the 
each solid curves. The corresponding standard model contribution is shown by the dashed 
line with its value shown along it. The Higgs mass is same as Figure 1. 




Figure 4: The tree level differential cross-section for e+e" ZHH at 1 TeV e+e" linear 
collider. The solid line correspond to ADD predictions for Ms — 2 TeV and n = 3, while 
the dashed line represent the standard model contribution, assuming Mh = 120 GeV. 




Figure 5: The fractional deviation R{ defined in tfie text) of the cross-section for e+e~ — > 
HHZ from the standard model prediction as a function of Ms for different choices of extra 
dimensions as shown along the curves. The solid and dashed curves correspond to y/s — 1 
TeV and 0.5 TeV NLC. The Higgs mass is same as Figure 1. 
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Figure 6: Variation of the cross-section (fb) for HH^ production with machine energy in the 
ADD model with Ms = 3.5, 4.5 and 5.5 TeV and n = 3. The Higgs mass is set to 120 GeV. 
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Figure 7: The tree level differential cross-section for e'^e — > HH'y at 1 TeV e'^e linear 
collider, M5 = 2 TeV, n = 3 and M^^ = 120 GeV. 
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Figure 8: Variation of the cross-section (fb) for HH^ production with Ms- The sohd, dashed 
and dotted curves with the extra dimension n = 3 — 6 from right to left correspond to ^/s — 3 
TeV, 1 TeV and 0.5 TeV respectively. 



